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Abstract 

The dispersion relation of energetic (few TeV) neutrinos traversing a 
medium is studied. We use the real time formalism of thermal field the- 
ory and we include the effects from the propagator of the W gauge boson. 
We consider then the MSW oscillations for cosmic neutrinos traversing the 
Earth, adopting for the neutrino parameters values suggested by the LSND 
results. It is found that the t'^ flux, for neutrinos passing through the center 
of the Earth, will appear reduced by 15% for energies around 10 TeV. 
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1 Introduction 



The quest to learn whether neutrinos are massive has been long and ardu- 
ous. In the recent years there are mounting astrophysical and laboratory 
data suggesting that neutrinos oscillate from one flavor to another,which can 
only happen if they have non-zero mass. The SuperKamiokande experiment 
on atmospheric neutrinos^], solar neutrino experiments^] and LSND accel- 
erator experiment [0] , provide three distinct scales of neutrino mass-squared 
differences: 



Amf,„ ~ 5 X 10-^eV^ 



Aml^ ~ 10-^eV^ (1) 



2 



^rnisNn - 10 eV 

To accomodate all data a mixing scheme with four massive neutrinos, 
three active (t'e, t'/^, z^r) and one sterile (z/^) is required. The mass pattern 
contains two pairs, each consisting of two nearly degenerate states, separated 
by the 'LSND gap' of a few eV|^. A natural assignment would have (f^, Ut) 
with a mass of few eV each and a mass-difference Am^^^ accounting for the 
atmospheric neutrino oscillations, while i^e and i^s, much lighter and with a 
splitting Am^^„, provide the MSW solution for the solar neutrinos. 

Few - eV mass neutrinos play a cosmological role as a hot dark matter 
(HDM) component in the popular cold-|-hot dark matter (CHDM) model0, 
providing less structure on small scales. In addition, neutrinos with eV 
masses mimic extra radiation at the epoch of matter-radiation equality and 
affect the cosmic microwave background (CMB) anisotropyP|. Sterile back- 
ground neutrinos may generate large neutrino asymmetries and influence the 
big bang nucleosynthesis (BBN)[0]. Massive relic neutrinos participate in 
gravitational clustering around galaxies. Recent calculations]^] indicate that 
gravitational clustering provides locally an increase in the neutrino density 
by a factor 10^, over the uniform density of the big bang cosmology. It has 
been suggested that the highest energy cosmic rays result from the an- 
nihilation of extremely high energy cosmic neutrinos on the background of 
gravitationally clustered relic neutrinos. From all the above it is clear that 
the issue of neutrinos with a mass of few eV is an interesting one and de- 
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serves further study. We suggest in this paper another approach to explore 
this issue. 

Solar neutrinos, in the MeV energy range, undergo resonant conversion 
with Aml^n ~ 10"^ eV^. We infer that the range ^m\g^^ ~ 10 eV^ 
can be studied more appropriately, via MSW oscillations, using neutrinos in 
the energy range of few TeV. The outer space provides powerful radiation 
sources, like Active Galactic Nuclei(AGN). AGN are the central regions of 
certain galaxies in which the emission of radiation can rival or even surpass 
the total power output of the entire galaxy by as much as a thousand fold. 
The source that powers AGN is believed to be gravity, i.e. matter accretion 
into a supermassive black hole located at the center of the galaxy. Accelerated 
protons may interact with matter or radiation in the AGN to produce pions 
whose decay products include neutrinos. It is anticipated that AGN could 
be the most luminous high energy neutrino sources in the universe and the 
diffuse isotropic neutrino flux from all AGN has been estimated |]10|. Also 
Gamma- Ray Burst (GRB) sources have been proposed as generators of high 
energy neutrinos |]10|, |lT|. 

Cosmic high energy muon neutrinos (E^, > 1 TeV) can be observed by 
neutrino telescopes |12|, |T3| by detecting the long range muons produced in 
charged current muon neutrino-nucleon interactions. The effective detector 
volume is enhanced in proportion to the range of the produced muon (typ- 
ically several kilometers). To reduce background, at the detection site one 
looks for upward moving muons induced by neutrinos traversing the Earth. 
We study then neutrino matter oscillations inside the Earth (z/^ ^ z/g) for 
neutrino energies above 1 TeV and with mass-squared difference and vacuum 
mixing angle for the neutrinos as suggested by the LSND results. 

The usual MSW analysis has previously been applied to low energy neu- 
trinos, where the Fermi model is an adequate description. Since we are 
interested in energetic neutrinos (Ej, > 1 TeV), the gauge-boson propagator 
effects might be important. In the next section we study the dispersion re- 
lation for energetic neutrinos in a medium, using the real-time formalism of 
finite-temperature field theory. In section 3 we apply our results to neutrino 
matter oscillations inside the Earth and present our conclusions. 
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2 Dispersion relation for neutrinos 



At finite density and temperature the properties of particles deviate from 
their vacuum values. This applies in particular to the dispersion relation 
which governs their plane- wave propagation. For neutrinos traversing a "fla- 
vor birefringent" medium, the effective mixing angle can become large, in- 
ducing a complete reversal of the flavor content of the neutrino state (MSW 
effect [ll4| ) . The real-time formalism of thermal field theory [|l^ is well suited 



to analyze such a problem. The real-time propagator consists of a sum of 
two parts - one corresponding to the propagator at zero temperature and the 
other corresponding to a temperature dependent part. The fermion propa- 
gator is given by 

Sfip) = {i> + + tV{p)] (2) 

r(p) = 27r5(/-m2)^(po)r;/(bo|) (3) 

where is the occupation number for the background fermions at temper- 
ature T = ^ and chemical potential /i. While the zero temperature part of 
the propagator can be thought of in the conventional manner as represent- 
ing the exchange of a virtual particle, the temperature dependent part, see 
eq.(3), represents an on-shell contribution. In a plasma of particles, there 
is a distribution of real particles and these real particles participate in the 
scattering process in addition to the exchange of virtual particles. 

Neutrinos propagating in a medium acquire an effective mass. Denoting 
by S the modified self-energy, the propagation of a neutrino in a medium is 
governed by the Dirac equation 

- = (5) 

At one- loop level S has the general form|p!6| 



E = m-(a^ + ^)i(l-75) (6) 

where m is vacuum neutrino mass, k is the neutrino four-momentum and u 
is the four-velocity of the medium. The coefficients a, h are functions of the 
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scalar quantities k'^, v? = 1 and k ■ u = u, with to the energy of the neutrino 
in the rest frame of the medium. The dispersion relation is given by 



det{ft - S) = (7) 
and to lowest order in small quantities provides 

- \k\'^ = rn^ - 2bco (8) 
Equivalently the refraction index n {\k\ — nuj) becomes 

n = 1 - — + - 9 

Neutrinos traversing the Earth interact with electrons via neutral current 
and charged current interactions. While all active neutrinos interact via the 
neutral current, only the z/g participate in the charged current interaction. 
The relevant one loop contribution to the self-energy of the neutrino is given 
by (see Fig. 1) 




Fig. 1 



.^^fj -^,ni-l,W + rn)iT{p)D,^{k-p) (10) 

where g is the SU(2) coupling constant and D^y is the propagator for the W 
gauge boson. In the unitary gauge we have 



^^i/(^ - P) 



—I 



{k-py-M^ 



{k-p)^{k-p)^ 



(11) 



Usually the energies involved are small compared to and the mo- 
mentum dependence in the W propagator is ignored. Since we are interested 
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in energetic neutrinos (E,^ > 1 TeV) we expand the propagator in powers of 



M^^. We obtain 

D^Ak - P) ^ + """''^""^"'''■'""'''1 (12) 

It is the above expression we have employed in our calculations. Notice 
also that from the electron propagator we kept only the term emanating from 
the medium. The first term in eq.(2), electron propagator in the vacuum, 
provides the wave-function renormalization of the neutrino in the vacuum (of 
no relevance to us). 

Working within the Fermi model approximation, i.e. keeping only the 
first term from the W propagator, eq.(12), we obtain for the corresponding 
self-energy Eo (evaluated in the rest frame of the medium) 

iE„ = iG^V2 7° \{l - 75) 2 / ^r;; (13) 

Since the number density of the background electrons is given by 

^. = 2/jgj,, (14) 

we arrive at the following expression for S^, valid in an arbitrary frame 

iE^^iCpV^iVe^ ^(1-75) (15) 

We read off then that 

6„ = -V2GFNe (16) 

in agreement with previous calculations [16-21]. 

Including the second term from eq.(12), i.e. the correction from the finite 
mass W propagator, we find for the corresponding self-energy Ei (evaluated 
also in the rest frame of the medium) 



iEi = —i 



• ^.7o^(i - 75) / + + ^ 



^^'"2^ ""J (27r) 



By terms we denote all these terms proportional to responsible for the 
a coefficient (see eqn.(6)). Since we are interested in finding the b coefficient 
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we ignore these terms, e and p are the energy and the momentum of the 
background fcrmions. Wc consider two distinct cases: relativistic electrons 
(R), where ~ and non-relativistic electrons (NR), where « e^. 
Then we have 

£ + ^~ce (18) 
with c— 1(1) for the R(NR) case. For Ei in an arbitrary frame we obtain 
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i9''_^.r /.I 



iEi = --f-^E,N,c <e> - 75) + ^ terms (19) 



with < e > the average energy of the electrons in the medium. We extract 
that 

h = 2v^cG'^7v /- > (20) 

Altogether 

b=bo + h^ -V2GpN,[l - 20 ^" ^ ] (21) 

With in the TeV range and electrons in a hot plasma the second term 
in the square bracket could be a significant correction. In our case the elec- 
trons inside the Earth can be considered as a Fermi gas and the maximum 
momentum pp is fixed by the local number density. For the actual Earth 
densities we find low values for p^. Therefore we obtain 

b^-V2G,N,[l-2^] (22) 

We deduce that the correction is not significant in the Earth and in our 
calculations we use the standard MSW formalism. 



3 Results and Conclusions 

Foccusing our attention into f ^ < — > oscillations, any neutrino state can 
be written as 

\u{t) >= ae{t)\ue > +a^{t)\l^^^ > (23) 
and the time evolution is given by 

4f"0=^("0 (24) 
dt \af^J \OLij,J 
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where 

A^]_( -Acos26'o + 2p7 Asin2eo \ 

Ap\ Asin2eo Acos29a - 2p-i ) ^ ' 

The neutrino momentum is denoted by p, is the vacuum mixing angle and 
A is the mass-squared difference. Matter effects are determined by 7||1^ 
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V2GfN^ (26) 



with Nf. the density of electrons in the medium. For the neutrino parameters 
we have chosen values implied by the LSND results 

A = 10 eV^ sin2eo = 0.07 (27) 

The relative population of cosmic neutrinos at the source is 2 : 1 u^. 
The smallness of the vacuum mixing angle suggests that vacuum oscillations 
(during the journey from the cosmic source to the Earth) do not alter the 
relative proportion. However during the passage through the Earth matter 
oscillations might change significantly the relative abundance of z/^ and Ue- 
The density of the Earth is accurately described by the two density model, 
where the core and the mantle each have a separate and constant density ||22 



Defining by (/) the angle between the neutrino direction and the tangent to the 
surface of the Earth, we considered neutrinos incident at = | (i.e. passing 
through the center of the Earth). In the absence of any oscillations, the 
muon neutrino flux exiting the Earth is identical to the initial incident muon 
neutrino flux. In the presence of oscillations, the muon neutrino flux to be 
detected, is equal to the initial incident muon neutrino flux times the factor 
[ P{Uf^ v^) + \P{h'e T^fj.)]- Fig. 2 presents this factor as a function of the 
energy. At around 10 TeV and for = ^, a 15% reduction in the number of 
muon neutrinos is observed. Notice that the observed dip corresponds to an 
MSW resonance. However at resonance the oscillation length is given by ||22 



^ 'ytan20o ^ ^ 

For small vacuum mixing angle 9o, becomes large and exceeds the radius 
of the Earth. For this reason the MSW dip is not fully developed. 

Neutrino telescopes will provide a new observation window on our cosmos 
and will help us probe the deepest reaches of distant astrophysical objects. 
The large energies involved might also allow us to enlarge our knowledge of 
particle physics. Determining the direction of the incident neutrino, we infer 
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the neutrino path through the Earth. Therefore the neutrino telescopes are 
suited for the study of neutrino oscillations in the GeV-TeV energy range, 
and for oscillation lengths of the order of the size of the Earth. Furthermore 
by comparing the event rates of neutrinos scratching the Earth, to neutrinos 
traversing the Earth, we obtain access to the neutrino oscillation probabil- 
ities. In our work we studied the oscillations among z/^ and Vf.- Another 
possibility is oscillations among and Pr and the appearance of Pr at TeV 
energies p3|, |2^ . This interesting case is under study. 
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